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Abstract 

Baryogenesis at the electroweak scale depends on divers but identifiable details of 
bubble wall dynamics and the particle physics. We show that inclusion of the dynamics 
of relative phase in two-Higgs-doublet model (2HDM) enhances the adiabatic order of 
the mechanism proposed by McLerran-Shaposhnikov-Turok-Voloshin where the scalar- 
scalar-vector triangle diagram with top quark in the loop gives rise to a significant 
contribution to the effective chemical potential biasing the Chern-Simons number. We 
also show that in 2HDM with less stringent constraints on Yukawa couplings than those 
imposed by natural flavour conservation, there are additional diagrams contributing to 
the effective chemical potential. These two effects can combine with several others to 
produce adequate baryon asymmetry at the electroweak scale. 

1 Introduction 



The main motivation for the models of CP violation has been to explain the e'/ e parameter 
for neutral kaon decay simultaneously preserving the vanishing smallness of the electric dipole 
moment of the neutron. The idea that CP violation can be accomodated in extensions of the 
Higgs sector was put forward by Weinberg in his three Higgs doublet model, and T.D. 
Lee [0] in his two Higgs doublet model (2HDM) with spontaneous CP violation. In these 
models the Yukawa couplings have to be constrained to suppress flavor changing neutral 
currents. A natural condition is to require that either (pi couples to the up quarks and 2 
to the down quarks or 0i couples to all the quarks and 4>2 does not couple at all. A discrete 
symmetry <f>\ — > —<\>\ is imposed on the Higgs potential to ensure suppression. In majority 
opinion the minimal standard model suffers from inadequate CP violation coming from the 
CKM matrix to be considered a viable model for baryogenesis at the electroweak scale. So 
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the 2HDM with with CP violation in the Higgs sector has attracted attention for the past 
several years in the context of electroweak baryogenesis. 

McLerran-Shaposhnikov-Turok and Voloshin had shown that in this model, the tri- 
angle diagram with one gauge boson leg and two scalar legs and top quark in the loop 
contributes a correction to the effective action ~ ON cs , where O is an operator that can be 
found out by explicit evaluation of the diagram and Ncs is the Chern-Simons (CS) number. 
This correction biases the CS number and since ANb = rif A Nqs, nett baryon number 
results, with the sign determined by the sign of < O >. 

Recently Yue-Liang-Wu f| has generalised the model of Glashow and Weinberg || and 
has proposed that the value of e'/e and the limits on EDM can be well explained with less 
stringent conditions on the Yukawa couplings. The Higgs boson-fermion couplings are now 
much more general and no discrete symmetries are imposed on the Higgs potential. In this 
class of models fermions acquire mass from both (pi and 02- Consequently in addition to the 
triangle diagrams with X on the external legs, diagrams with 2 also on the external legs 
start contributing to the effective action. 

The MSTV proposal has been criticised on the grounds that it is a higher order adiabatic 
effect, of the order ((0) T /T) 4 . In their considerations the relative phase between the two 
Higgs vacuum expectation values was treated constant throughout bubble evolution. Here 
we provide the details of the argument || that the inclusion of the dynamics of the relative 
phase removes the adiabatic suppression, making it an effect of the order ((</>) T /T) 2 . We have 
also previously shown [0] that the uncertainties of physical parameters of bubble formation 
H can be circumvented if we consider a string induced phase transition 0). This is the 
scenario we shall be considering here as well. It may be noted that since according to 0, the 
string induced bubbles provide adiabatic conditions, all B-genesis mechanisms relying on such 
conditions are workable. In particular the mechanisms of Cohen and Kaplan |T(J and Cohen- 
Kaplan-Nelson jri] can also proceed through string induced bubbles. We recapitulate here 
the corresponding bubble profile ansatz including the time evolution of the relative phase and 
calculate the resultant contribution to the effective action. Putting together the two sources 
of enhancement, viz., dynamics of the relative phase and additional triangle diagrams, we 
show that the mechanism has sufficient potential to give adequate baryogenesis. 



2 The triangle diagram and correction to the effective 
action 

It was first pointed out by Turok and Zadrony [T!| and then explicitly calculated by 
McLerran-Shaposhnikov-Turok and Voloshin ||, that for 2HDM a term biasing the Chern- 
Simons number with a CP odd chemical potential is contributed by the triangle diagram of 
figure 1. This contribution to the effctive action at finite temperature is 

AS = ^ C (3)(^) 2 T |^^ x J (V^Vofa+V^V^e^F^x (1) 
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where m t is the mass of the top quark, ( is the Riemann Zeta function, and a a are the Pauli 
matrices. For homogeneous but time varying configurations of the Higgs fields and in the 
A a = gauge, 

ASx = ^C(3)(^) 2 ^/ dttfJVofa - Voti^Ncs (2) 

= O^cs (3) 

Clearly, here the top quark acquires mass only from 0i as demanded by the Glashow- 
Weinberg natural flavor conservation (NFC) criteria. In Wu's model, the FC criterion is 
satisfied under relaxed conditions on the Higgs-fermions couplings. In this case the general 
Yukawa interactions can be written as 

Ly = qiL{T a D ) tj D jR <pa + qlLiTuhjUMa + kL{T a E ) tj E ]R <p a + h.c (4) 

where k and0 a are SU(2)l doublet quarks, leptons and Higgs bosons, Ui, D iy Ei are 
SU(2)l singlets. % — 1, 2, ...n is a generation label and a = 1, 2 is a Higgs doublet label. TV 
(F = U,D,E) are the Yukawa coupling matrices. According to Wu the Glashow- Weinberg 
criteria can be replaced by a theorem which states that the flavor conservation for the neutral 
currents is natural in the Higgs sector or equivalently, the matrices T a F (F = U,D,E) are 
diagonalizable simultaneously by a biunitary or biorthogonal transformation, if and only if 
the square n x n T a F are represented in terms of the linear combinations of a complete set 
of n x n matrices {Vt a F) a = 1, 2, ...n). 

r a F = T,9 F aa^F a (5) 

a 

where, £lp a satisfy the following orthogonal condition f2^ a (f2^ /3 )^ = L a p5 a i3, 
(Qp^flpP = R a p5 a/ 3 , with the normalization J2 a L a F = H a R a F = 1- This generalised 
Yukawa coupling prompts us to consider another set of triangle diagrams as shown in fig- 
2. The contribution at finite temperature from these diagrams to the effective action can 
similarly be calculated to be 

AS 2 = Z!Zc(3)(-L) 2 r 22 | dt[^{V^ 2 ) - (Vo^^Ncs (6) 

AS 3 = ^C(3)(^) 2 r 1 r 2 y dt[^(V <h) - (VohYHNcs (7) 
where T^) is the Yukawa coupling when the fermions couple to 0i(2)- 



2.1 The bubble profile 

Now to find the bubble profile we use the following ansatze for the finite temperature vacuum 
expectation value of the Higgs fields 

0!° = pi(r,f)e-* (t) (8) 
02° = P2(r,t)e^ (9) 
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where, as pointed out by Cohen-Kaplan-Nelson [O we can fix the unitary gauge ensuring 
9 is the physical pseudo-scalar orthogonal to the Goldstone boson eaten by Z. This gauge 
fixing gives the relation between oj and 9 to be, 

d„u = (pi/p 2 ) 2 d ll 9 (10) 

We assume that the phase transition takes place when a combination of p\ and P2 becomes 
massless in a particular direction 7 in the <\>\ — 02° plane. Hence, in the bubble profile we 
may take 

pi(r,t) = p(r,t) COS7 P2(r,t) = p(r,t) sin 7 (11) 

with this parametrization and finite temperature corrections, the effective potential for the 
p, 9 and u is 

V T (p,9,u) = M x \T)p 2 -ETp* + ^p' 

+M 2 2 (T) (cos i cos (u + 9) + Ci sin f sin (0 + cu))p 2 
+yp 4 (cos 2 (# + ^) + C'ism 2 (# + ^)) (12) 

where, M\ 2 and M2 2 are temperature corrected mass parameters and K\ and K 2 are com- 
binations of quartic coupling constants with small temperature dependent corrections. The 



constant C\ is the ratio of A5 and A6 in the standard parameterisation of the 2HDM [[14 
and £ is the phase of the m 
of (j>2- Rescaling p — > (2E' 
potential can be written as 



and £ is the phase of the neutral component of zero temperature vacuum expectation value 

™ and t -> tr-W 



of 2 . Rescaling p — ► (2ET/K 1 )p / , r — ► r' ' j^W- and t — ► t' /^W- omitting the primes the 



Vr(p,8,u) = (^) 4 [^(l + Bl )p 2 -f/> 3 + ^ 4 
K 

+ £ 2 )(cos£cos(cj + 9) + Cisin£sin(6> + cj))p 2 
+^p 4 (cos 2 (fl + u) + d sin 2 (# + w))] (13) 

where E 1 = — 1 an d -E2 = ~~ ^" ^ e g eome t r y of string induced bubbles is 

cylindrical. With this in mind, the time dependence of p and 9 can be found by solving the 
following equations, 

d 2 p d 2 p ldp (d9\ 2 3 r ^ 2 ^ 3 K lM _ 

+K 2 d sin 2 (9 + Op 3 + -y (1 + E 2 ) (cos B + C x sin £ sin (f + B))p = (14) 

^ + ^sin2(^ + 0p 2 

+^(1 + E 2 )(C sine cos (0 + -sin(0 + o;-O) =0 (15) 



The time independent solution of p can be found by imposing the boundary condition p — > 
as r — > oo and % — ► as r — > 0. Subsequently, as the nontrivial minimum becomes 
favorable, the same solutions begin to evolve in time. The parameters used in the equations 
are, 

Kx, K 2 = 0.1 - 0.001, 

d = 1.1, C = C x -l,Ex = -0.074, E 2 = -0.07, 
£ = 0.2, E ~ 0.01 

The time evolution of the bubble profile p(r,t) has been reported earlier 0. In fig-3, we 
show the time evolution of the relative phase in the regions where p has become nonzero. 
If the initial reference value is zero, it oscillates to reach the stationary value dictated by 
the 2HDM effective potential at the relevant temperature. We assume this value to be 0(1) 
since no natural reasons prevent it from being so. 

2.2 Evaluation of the operator 

Now we can use these solutions to evaluate the average value of the operators as, 

1 + 2 = 28C(3) (-J j^Ai J p 2 ^dt (16) 

and, 

Q 3 = 28C(3)Q 2 fU/ P'f/^dt (17) 

where A\ = (PiCOS7) 2 + (r 2 sin7) 2 and A 2 = TiT 2 sin2 / y. These parameters have to be 
determined from the phenomenology of the 2HDM which is yet far from being tested by 
current experiments. Knowing the solutions to eqn.s (14)-(15), we can estimate the integral 
in eqn. (16) to be p 00 2 A9 where A9 is the nett change in the relative phase 9 at any given 
point as the bubble wall sweeps past it, and = 1. There is an additional contribution 
from the transient part of 9, which can also be calculated numerically, but is not significant. 
As for the term (9 3 , its CP odd part has the magnitude / o 2 XD A(sin6 l ) from arguments already 
given. Putting in other known factors, we see that 

O l + 2 ~ A 1 (E/K 1 fA9 (18) 

and a similar contribution from (9 3 . Recall that E is the dimensionless cubic self-coupling 
induced by thermal loops, and K\ involves the quartic self-couplings of the 2HDM. For 
naturalness we would like A9 to be 0(1) but the remaining factor is numerically a small 
magnitude, perhaps between 1 and 10~ 4 . Note however that the effect is not suppressed by 
the physics of the process viz., the bubble wall dynamics. We emphasise again that this is 
the consequence of the dynamics of the relative phase of the 2HDM. 
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3 Estimation of the asymmetry 



To estimate the baryon asymmetry we assume the presence of high temperature sphaleron 
processes inside the bubble wall. The rate of such transitions per unit time per unit volume 
is of the order ~ Ka w T 4 , k ~ 1 |R) |. The number of fermions created per unit time in the 
bubble wall is given by 

B = K(a w T) 4 lS x (19) 

where we have made use of a well established master formula ]16| [|13|], and where / and S 
are the thickness and the surface area of the bubble wall respectively. From which we get 
the baryon to photon ratio to be 



n B 1 



u: ) 



Neff 

E ,2 



1(T 8 x ( j^- ) A9 (20) 
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where we have used a w ~ 10~2 and N e ff ~ 100. This answer easily accomodates the 
observed value of this number. 

It is worth emphasising the physics of this answer which is fairly robust against changes 
in the specific particle physics models. The thermal rate contributes 10 -6 through a^, and 
another 10~ 2 is contributed by N e ff. The remaining smallness of the answer follows from 
smallness of the thermal vacuum value (<ft) T ~ ET/K 1 , which is small on the scale set by 
T c . The appearance of this particular physical quantity has to do with our picture of the 
process occuring in bubble walls as the phase transition is yet in progress. 



4 Conclusion 



The MSTV proposal has been criticised jn| on the grounds that the operator O is of the 



order of (0/T) 4 . Since <fi, i.e., the temperature dependent vacuum value (0) T is T, the 
effect was thought to be unacceptably small. But with time variation of the relative phase 
allowed, we find the relevent operator is of the order of (0/T) 2 . The effect is therefore not 
intrinsically suppressed. We also note that the same analysis could be fruitfully applied to 



the proposals of |I(J and |fLl]j . Secondly, we have shown that in principle there are more 
diagrams in the 2HDM contributing to the effective action. Although in making numerical 
estimates we are hampered by large number of unknown parameters in the 2HDM, it is 
worth remembering the existence of these effects as potential sources of enhancement in 
electroweak baryogenesis. 
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Figure captions 

• Fig-l The contribution to S e ^ from one higgs. 

• Fig-2 Additional diagrams contributing to S e ^ for generalised yukawa couplings. 

• Fig-3 Evolution of the relative phase during passage of the bubble wall. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9508330vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9508330vl 



